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Abstract. We consider filamentary and planar large scale 
structures as possible refraction channels for electromag- 
netic radiation coming from cosmological structures. By 
this hypothesis, it is possible to explain the quasar lu- 
minosity distribution and, in particular, the presence of 
"twin" and "brother" objects. Method and details of sim- 
ulation are given. 



1. Introduction 

Gravitational lensing (Straumann et al. 1998) is today 
. considered a basic tool to investigate the large scale struc- 
' ture of the universe and to perform tests of cosmological 
models. Furthermore, it plays a relevant role in the dark 
matter search from small to large scales. On the other 
hand, the lensing effects can be explained by considering 
the action of a weak gravitational field deflecting the light 
, rays which come from a fixed source. In other words, it 
is possible to define a refraction index, connected to the 
, Newtonian gravitational field. 

' However, it is worthwhile to note that, in optics, there 
are other kinds of effects and instruments, like optical fi- 
bres and waveguides, which use the same deflection prin- 
ciples: The analogy with the gravitational field can be ex- 
tended also to these processes taking into account a sort 
of " nonlinearity" in gravitational lensing. 

Considering the concept of "channel" of sister sciences, 
we have refraction channels, acoustic channels and several 
other examples where propagating waves are trapped by 
some structures. 

In cosmology, matter distributions, able to produce 
gravitational potential inducing such "channelling ef- 
fects" , could explain several anomalies which are hard to 
connect to standard gravitational lensing models, like, for 
example, the existence of objects with the same spectrum 
and redshift, but placed in the sky at large angular dis- 
tances (e.g. twins). The huge luminosity of quasars (Fan 
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et al. 2000; Zheng et al. 2000) could be explained by using 
planar cosmological structures too. In fact, taking into 
account expanding bubbles, produced during first order 
phase transition, a matter density gradient is generated 
when they come in touch each other. Such expanding bub- 
bles create tridimensional spatial structures like hexahe- 
drons (e.g. honeycombs) to minimize the interstitial space. 
However the growth is different with respect to crystal 
models since in these last cases the growth is stopped in 
such a direction as well as a cell of the lattice touches with 
an other expanding crystal in the same direction. In cos- 
mology, we have to take into account the expansion of the 
universe. 

Our hypothesis is that the surfaces of such hexahe- 
drons, which are contact surfaces, could be approximated 
by planar structures acting as channels for electromagnetic 
radiation. 

Furthermore, as we will discuss below, also filamentary 
structures could act like channels too. 

In principle, all primordial topological defects as 
cosmic strings, textures and domain walls (Vilenkin 
1984; Vilenkin 1985) could be evolved so that today 
they can trap electromagnetic or gravitational radiation 
(Capozziello & lovane 1998; Turner et al. 1986). 

This gravitational waveguiding effect has the same 
physical reason of the gravitational lensing effect, which 
account for the light deflection given by a gravitational 
field acting like a refraction medium . However, there are 
some substantial differences: the gravitational lenses are 
usually compact objects. On the contrary, the density can 
be low in these gravitational channels; moreover the radi- 
ation is deflected by the lens, but it does not pass through 
the object, while the light have to travel "inside" the mat- 
ter distribution in the case of gravitational channels. 

The existence of a waveguiding effect could explain the 
huge luminosity of very far objects like quasars. Another 
effect could be the generation of twins and brothers. In 
other words, an observer could detect two images of the 
same object: The real image and the one catched by grav- 
itational channel. 
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From a theoretical point of view, the matter distri- 
bution should locally produce an effective gravitational 
potential of the form $ (r) ~ r^. 

It is worthwhile to note that the specific form of r (and 
then of the effective gravitational potential) depends on 
the particular symmetry of the situation. As we shall dis- 
cuss below, wc have r = x for " planar" mass distributions 
while r = \/3p+y^ in the case of " filamentary" distribu- 
tions. In both cases, we are assuming an almost constant 
density p which is substantially different with respect to 
the background. 

The planar or filamentary distributions, essentially 
connected to a gradient of matter density, could be gen- 
erated at the contact surfaces among several expanding 
bubbles. 

The paper is organized as follows. In the Sec. 2, the 
theoretical model for gravitational channels is developed. 
Sec. 3 is devoted to the discussion of early cosmological 
phase transition, topological defects and possible connec- 
tions to the today observed large scale structures. In Sec. 
4, we present a simulation by which we try to explain the 
luminosity distribution of quasars assuming that a part 
of them has been " twined" because the radiation emitted 
from them have undergone a channelling process by inter- 
vening filamentary or planar large scale structures. The 
results are discussed in Sec. 5. Conclusions are drawn in 
Sec.6. 



Mr) ( $(r)\ 

500^1 + 2^; ~ -^ife h - 2^ j , (3) 

where, $(r) is the Newtonian potential and we are 
assuming the weak field, 1, and the slow motion 

approximation \v\ <C c. 

Channeling solutions can be obtained, by reducing 
Maxwell equations in a medium to the scalar Helmoltz 
equation for the fields E and H . For the electrical field, 
after Fourier transformation, we can write (Capozziello & 
lovane 1998) 



AE,,(r) ^n^(w,r)E^(r) = 



(4) 



where io is the frequency of the given Fourier compo- 
nent. 

A similar equation holds for the magnetic field. The 
general solution has the form 



E^(r)=E£0)(r)+ / e(f,r',a;)V 



ES')(r)Vg(u;,r) 
£(a;,r) 

where Q is the Green function for the equation 



(5) 



2. Light propagation by gravitational refraction 
channels 

The behaviour of the electromagnetic field, without source 
and in the presence of a gravitational field (Schneider, 
Ehlers, Falco 1992) can be described by the Maxwell equa- 
tions 



dx-i 



+ 



dF, 



Pi 



dx° 



+ 



dF. 



7Q 



dx0 



0; 



1 d_ 



gF"^) = 0, 
(1) 



where F"'^ is the ckictroniagnetic field tensor and ^/—g 
is the determinant of the four-dimensional metric tensor. 
For a static gravitational field, these equations can be re- 
duced to the usual Maxwell equations describing the elec- 
tromagnetic field in media where the dielectric and mag- 
netic tensor permeabilities are connected to the metric 
tensor g^,^ by the equation 



£ifc = Mife = -W^^[det5fife] ^^^Qik i,k= 1,2,3. (2) 

If the medium is isotropic, the metric tensor is diago- 
nal and the refraction index is n (r) = {efi) ^ (Capozziello 
& lovane 1999). 

In the weak field approximation, the metric tensor 
components are (Blioh & Minakov 1989; Landau & Lif- 
sits 1975) 



a(r,r',a;) = 5(r-r'). 



(6) 



This is a Schrodinger-like equation for the energy con- 
stant E = 0. In fact, if we write down the Hamilto- 
nian operator Ti. = — -|- J7(r) ,with h = m = 1 and 



1^2 

U(r) = — 2— n^(w,r) . a fully analogv between the two 

problems is obtained. U{r) assumes the role of the New- 
tonian potential ^'(r). 

Let us consider the scalar Helmoltz equation for an 
arbitrary monochromatic component of the electric field 



d^E d'^E d^E ,22,, 



(7) 



where k is the wavenumber. The coordinate z can be 
considered as the longitudinal one and it can measure the 
space distance along the structure produced by a mass 
distribution with an optical axis. 

Let us consider now a solution of the form 



E = Hq ^^^^ exp J no{z')dz'^ ; no = n(0, 0, z) , 

(8) 

where ^'(a;, y, z) is a slowly varying spatial amplitude 
along the z axis, and exp(ifcnz) is a rapidly oscillating 
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phase factor. Its clear that the beam propagation is along 
the z axis. We rewrite the Helmoltz equation neglecting 
the second order derivative in longitudinal coordinate z 
and obtain a Schrodinger-like equation for 



- 



(9) 



where A is the electromagnetic radiation wavelength 

dz' 

and we adopt the new variable ^ = I — is a new 



no{z') 

normalized variable with respect to the refraction index. 
It is interesting to stress that the role of Planck constant 
is assumed by A. 

At this point, it is worthwhile to note that if one 
has the distribution of matter in the form of a cylinder 
or a sphere with a constant (dust) density, the gravita- 
tional potential inside has a parabolic profile providing 
channelling effects for electromagnetic radiation analo- 
gous to the sel-foc optical waveguides realized for opti- 
cal fibres. In this case, Schrodinger-like equation is a two- 
dimensional quantum harmonic oscillator which solutions 
provs the form of Gauss-Hermite polynomials (see, for ex- 
ample (Man'ko 1986)). 

The channeling effect depends explicitly on the shape 
of the potential: The radiation from a remote source does 
not attenuate if t/ ~ ; this situation is realized suppos- 
ing a "planar" mass distribution (r = x) or a filamentary 
distribution (r = \/3p+y^) with constant density . In 
other words, if the radiation, travelling from some source, 
undergoes a channelling effect, it does not attenuate like 
as usual, but it is, in some sense conserved; this 
fact means that the source brightness will turn out to be 
much stronger than the brightness of analogous objects 
located at the same distance (i.e. at the same redshift Z) 
and the apparent energy released by the source will be 
anomalously large. 

To fix the ideas, let us estimate how the electric field 
propagates into an ideal filament whose internal potential 
is 



AnGp 



(10) 



A spherical wave from a source has the form 



E={l/R)e^^{ikR). (11) 
In the paraxial approximation, Eq. (11) becomes 



E{z,r) 



1 



exp I ikz + 



ikr 



r 

2^ 



(12) 



where we are using the expansion 



R={z'' + r^) 



1 + 



2^2 



(13) 



It is realistic to assume no — 1 so that (, = z. 

Let us consider now that the starting point of the grav- 
itational channel of length L is at a distance I from a 
source shifted by a distance a from the structure longitu- 
dinal axis in the x direction. The amplitude of the field E, 
entering the channel, is 



1 

7 exp 



ikl-l 



y 



(14) 



,1/2 



so that we have R = (^P + y"^ + {x — a)^)"'" . We can 
calculate the amplitude of the field at the exit of the fila- 
ment by the equation 

*o«t {x, y,l + L) = 

= J dxidyiGix, y,l + L, si, yi, Z)*i„(a;i, yi, 0,(15) 

where G is again the Green function for ^. In particu- 
lar, taking into account the potential (10), G is the prop- 
agator of the harmonic oscillator; then "if out is a Gaussian 
integral that can be exactly evaluated. There are two in- 
teresting limit for "^out- 

If w/ <C 1, we have 



1 



exp 



which means that the radiation emerging from a point 
with coordinate (a, 0, 0) is focused near a point with coor- 
dinates (—a, 0,1 + L) (that is the radius has to be of the 
order of the wavelength). This means that, when the beam 
from an extended source is focused inside a gravitational 
channel at a distance L, an inverted image of the source is 
formed, having the very same geometrical dimensions of 
the source. The channel "draws" the source closer to the 
observer since, if the true distance of the observer from 
the source is R, its image brightness will correspond to 
that of a similar source at the closer distance 



R^ff = R-l-L. 



(17) 



In the opposite limit, loI ^ 1, we have tanujL oo, 
so that L ~ n/u!, that is the shortest focal length of the 
gravitational channel is 



^foc 



(18) 



It is relevant to note that the expression L/oc leads 
to values of about 100 Mpc for a typical galactic density 
Po ^ 10^^"' g/cm?. In other words, this can be assumed as 
a typical length for large scale structures. 
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3. Are there candidates for cosmological 
refraction channels? 

Early structures, resulting from inflationary phase transi- 
tions, could be good candidates to implement the above 
filamentary or planar structures acting as cosmological re- 
fraction channels. Essentially, if the primordial universe 
undergoes a first or second order phase transition, such a 
transition takes place on a short time scale ( t < H^^ ), 
and it will lead to correlation regions inside of which the 
value of a certain scalar field (the inflaton which leads 
dynamics) is approximately constant, but outside Lp ranges 
randomly over the vacuum manifold. In the simplest case, 
for a Ginzburg-Landau-like potential of the form 

V{^) = \{^^-^1)\ (19) 

the correlation regions are separated by domain walls, 
regions in space where leaves the vacuum manifold M 
and where, therefore, potential energy is localized. This 
energy density can act as a seed for structure. 

There arc various types of local and global topological 
defects (Kibble 1976) (regions of trapped energy density) 
depending on the number of real components of ip (being, 
in general, = 'Ym=i "^f)- ■'■^ domain walls 

n = 1. 

The rigorous mathematical definition refers to the ho- 

motopy of M . The words "local" and "global" refer to 
whether the symmetry, which is broken, is a gauge or 
global symmetry. In the case of local symmetries, the topo- 
logical defects have a well defined core outside of which ip 
contains no energy density in spite of non vanishing gra- 
dients V<^ ; the gauge field can absorb the gradient, i.e. 
D^if = when d^iip ^ 0. Global topological defects, how- 
ever, have long range density fields and forces. 

Practically, the topological defects can be classified by 
their topological characteristic which, in term of the scalar 
field if , depends on the number of real components. 

Pointlike defects (n = 3) give monopoles, linear defects 
(n = 2) give cosmic strings, smfaccs-likc defects {n = 1) 
are domain walls, hypersurface-like defects, but not prop- 
erly topological defects, with n = 4, are textures. 

In principle, all these structures can be seeds for large 
scale structures depending on the cosmological model we 
adopt. Taking into account realistic scenarios, the best 
candidates seem to be cosmic strings. 

Furthermore, any first-order phase transition proceeds 
via bubble nucleation so that invoking filamentary and 
planar structures for channeling effects is reasonable (local 
surface portions of large bubbles can be approximated to 
planar structures). 

For example, cosmic strings easily realize the condition 
to obtain a Newtonian potential useful, as discussed above, 
to produce a gravitational refraction channel. In fact, in 
the weak energy limit approximation, these topological 
defects are described internally by the Poisson equation 



= Pq and externally by = , where po is a con- 
stant. The condition on the density gradient is naturally 
recovered. Another useful feature is that they act as grav- 
itational lenses after the quasar formation (Vilcnkin 1985; 
Gott III 1985) Immediately we get, internally, $ ~ and, 
by the dynamical evolution of the universe, strings can 
evolve into structures with lengths of the order ~ 100 
Mpc. However, also after evolution, strings remain "wires" 
and then, in order to get the cylindrical structures which 
we need for gravitational channels (e.g. a filament of galax- 
ies) , we have to invoke further processes where strings are 
seeds to cluster matter at large scales. 

Furthermore, the motion of cosmic strings with respect 
to the background produces wakes and filaments which 
are able to evolve into large scale structures like clusters 
or filaments of galaxies (Vachaspati & Vilenkin 1991). 

All these arguments, as it is well known, are hypotheses 
for large scale structure formation. 

Our issue is now to construct reasonable toy models 
which, by channeling effects, could explain, in a simple 
way, the huge luminosity of quasars and twin effects in 
gravitational lensing. 

4. The model 

The major aim of our simulation is to perform a test about 
the possibility that twin or brother quasars could be ex- 
plained by trapped and guided light, which was emitted 
by quasar sources and guided thanks to structures acting 
like channels. Moreover, the single quasar itself could be 
a typical object (a bulge of a standard galaxy) which ap- 
pears so magnified thanks to the presence of channelling 
effects. 

In our simulation, as discussed above, we can take into 
account a gravitational first order phase transition from 
which false vacuum bubbles nucleate inside a true vac- 
uum cosmological background. Our fundamental assump- 
tion is that they cosmologically evolve giving rise to a cel- 
lular structure capable of yielding the observed large scale 
structures (Kolb & Turner 1990). After, we take into ac- 
count also filamentary large scale structures. 

It is also important to note that, due to the cosmo- 
logical expansion, the bubbles continue to grow even after 
they collide; this is the difference with models used in 
crystal growth such as the Voronoi and the Johnson-Mehl 
models, where crystals stop to grow in the direction of the 
neighbour one. 

By simple geometric arguments, we can put around a 
circle other six circles in a given area, with the same ra- 
dius. If these circles are expanding, to minimize the space 
among them, they degenerate into hexagons. The Fig. 1 
shows the situation. 

In other words, density gradients of matter will be 
generated among the hexagons in the boundary regions. 
In three dimension, the spheres will degenerate into do- 
decahedrons by forming a typical honeycomb structure 



S. Capozziello & G. lovane: Quasar luminosity and twin efTects induced by filamentary and planar structures 



5 




Fig. 1. Model for the formation of contacting planar struc- 
tures. 




Fig. 2. Elementary cell like dodecaedron. 

with 2-dimensional elementary pentagonal cells (see Fig. 
2), whose surfaces will give rise to planar structures with 
matter density different from zero. 

Finally, given an expanding universe several dodec- 
ahedrons, remnants of a primordial gravitational phase 
transition, could yield a sort of honeycomb structure. The 
isotropy and homogeneity would be preserved at larger 
scales (size ^ 1000 Mpc). 

In this case, the formation of clustered galactic struc- 
tures will be naturally due to the contact among several 
spheres, which would give rise to density gradients. 

Such planar structures on contact surfaces could act 
as refraction channels, e.g. they could trap and guide the 




Fig. 3. Elementary cell like trapezoidal icositetrahedron. 

light due to the quasi-constant surface density; in this way, 
we should see a universe like that observed. To be more 
precise, we do not see an exact honeycomb, but a selfsimi- 
lar structure (Sylos Labini, Mountuori & Pietronero 1998) 
like those described in fractal models: in a more realistic 
model, we have to consider several spheres with different 
radii Ri (we have not a perfect dodecahedron, but a more 
complicate polyhedron, see for example. Fig. 3). Dynam- 
ics will be governed by the symmetry breaking during the 
phase transition: vacuum density could be different in dif- 
ferent bubbles. However, we can obtain relevant results 
just taking into account balls with the same radius, that 
is, by considering bubble nucleation at same density and 
expanding rate at constant speed v. Similar considerations 
can be done also for filamentary structures which act as 
refraction channels. 

4-1- Hypotheses and previsions of the simulation 

Let lis consider a universe model with an observed radius 
i?o = 3000 Mpc. Furthermore, let us assume 10'* quasars 
uniformly and isotropically distributed in this spherical 
space. The positions of the objects are fixed in a random 
way by using a random number generator. Moreover, the 
quasars are point-like objects. Furthermore, it is supposed 
that every single quasar emits light in an isotropic way. 

Between the refraction channels, we have to distinguish 
filamentary structures and planar structures. 

We produce randomly distributed filamentary struc- 
tures with a dimension of about 100 Mpc along the major 
axis and a shape like a cylinder with 5 points, where the 
longitudinal axis of the cylinder changes its direction. A 
typical simulated channel is shown in the Fig. 4. 
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In particular, the channels of the model are built by 
five cylinders with a length of 1 20 Mpc and radius 
r ~ 1 100 Kpc. The positions of the channel starting 
points, that is the centres of the first base circle in the 
first cylinder and the longitudinal axis of the cylinder are 
randomly generated. 

There is a constraint on this generation: we assume 
that the cylinders cannot collapse in a compact object. In 
other words, by defining 6i and the angles in the links 
between two consecutive cylinders along their longitudinal 
axis (see Fig. 5), we impose: 

< -di < 7r/2 and < (fi < 27r, 

with -di and ipi randomly selected. 

Alternatively, we can have casually distributed planar 
structures acting like refraction channels, with linear di- 
mensions of about 100 Mpc, on pentagonal surfaces. 

In particular, in a universe with a radius Rq = 3000 
Mpc, we can inscribe several pseudo-balls bi with a radius 
Ri such that the linear dimensions have an extension of 
about li = 100 Mpc. This can be made in different ways: 

1) Let us consider a circular sector with an angular 
spread of 7r/24 = 7.5°; by assuming that the linear exten- 
sion of the structure is of 100 Mpc, we obtain: 

100 Mpc = li = Ri = 764 Mpc. (20) 

In this case, we can inscribe ntaiis = 60 balls with a 
radius Ri in the universe with a radius Rq. Moreover, we 
have to calculate the error by assuming 

27ri?i = 4800 Mpc, I = 100 Mpc ^ e = = 2%. 

zwRi 

2) By assuming a circular sector with an angular 
spread of n/18 = 10°, we obtain 

Ri = 573 Mpc, nbaiis = 143, e = 3%. 

3) By assuming a circular sector with an angular 
spread of 7r/12 = 15°, we obtain 

Ri = 382 Mpc, nbaiis = 484, e = 4%. 

4) By assuming a circular sector with an angular 
spread of tt/IO = 18°, we obtain 

Ri = 318 Mpc, nhaiis = 840, e = 5%. 

5) Finally assuming a circular sector with an angular 
spread of 7r/9 = 20°, we obtain 



Ri = 286 Mpc, ribaiis = 1154, e = 6%. 

In the simulation, it is useful to implement an attenu- 
ation mechanism to make the model more realistic. If the 
sources are at distances of the order of 1000 Mpc, since the 
luminosity decrc^ascs as R^^, we assume that only fluxes 
large enough to survive to this physical cut can be de- 
tected as the channelling effects intervene. Otherwise they 
are completely absorbed by intergalactic media. 

Next step is to verify how much light is trapped by the 
intervening channels and how many sources are lensed and 
then appear as twin images to the observer. 

4-. 2. The implementation 

The simulation is made by using the New Object Visual 
Programming Technique (National Instruments: Reference 
Manual 1998; National Instruments: Functions Manual 
1998; National Instruments: G-Math Manual 1998). 

The simulation is built up by a main program, that 
calls and manages a collection of subprocesses described 
as follows. 

1) The first process is built up by 5 secondary se- 
quences: 

1.1) The first one - sequence 1.1 - gives the quasar 
coordinates which are randomly generate; 

1.2) The second one - sequence 1.2 - gives the 
photons coming out from quasars: in this step the arrival 
points of coming out photons are evaluated; 

1.3) The third one - sequence 1.3 - generates the 
light rays at each quasar position; 

1.4) The fourth one - sequence 1.4 - rejects all 
photons out of 2i?o = 6000 Mpc (the boundary of our toy 
universe); 

1.5) The fifth one - sequence 1.5 - prepares a mul- 
tidimensional array and a file to perform and to save the 
intersections between the rays and the waveguides; 

2) The second process regards the objects, whose light 
is picked up, trapped and guided by the structures (it gives 
the twined quasars); 

3) The third process has two sequences: 

3.1) generates the refraction channel in random 
positions into the spherical volume representing the uni- 
verse; 

3.2) evaluates the parameters to intersect the rays 
from the quasars (this subprocess is completely equivalent 
to 1.5); 

4) The fourth process evaluates the attenuation func- 
tion for the light coming from quasars; 

5) The fifth process prepares the trapping; 

6) The sixth process performs the trapping and guiding 
of light in the channels; 

7) The seventh process tests that the same light rays 
is not counted more than one time, if it intersects more 
sub-cylinders of the same channel. 
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IppMpc 



Fig. 4. A typical simulated channel. 



200 Kpc 




Fig. 5. Scheme of a channel. 



In Fig. 6, there are shown several examples of simu- 
lated refraction channels. 



5. The results 

The results of the simulation are summarized, by consid- 
ering: 1) a fixed number of channels chosen in relation to 
Ri for the planar case and 200 for the filamentary one; 2) 
10'* quasars randomly fixed as input. 

a) The quasar positions R (in Mpc) and the redshifts 
after the attenuation, have a Gaussian distribution 

fit with the mean value and the standard deviation in 
full agreement with theoretical and observative results 
(Capozziello & lovane 1999; Fort & Mellier 1994; Broad- 
hurst, Taylor, Peakock 1995); 

b) Also for the twined quasars, R and Z are compara- 
ble with the theoretical models; in particular, the twined 
quasars are a subset of the input quasars and they are also 
isotropically distributed in the planes x — y and x — z. 

c) The Table 1 summarizes the obtained results for the 
planar channels from a numerical point of view. 

About 10% of initial quasars pass trough the atten- 
uation process. A fraction, from 7% to 35% of quasars 
after the attenuation, undergoes a channelling effect. This 



fraction is linked to the number of channels. The num- 
ber of twined quasars increases as the number of channels 
increases, by fixing the number of quasars after the atten- 
uation. A preliminary fit can be made between the num- 
ber of refraction channels and the twined quasars. Fig. 10 
shows a linear correlation. 

The results are in full agreement with the filamentary 

channels. In fact, in this case, we have 982 quasars after 
the attenuation and, among them, 137 are twined. In other 
words, about 10% of the initial quasars pass trough the 
attenuation process and a fraction of about 14% of these 
ones undergoes a channeling effect. 

In a first approximation, a superposition of the two 
distributions ( not twined quasars after attenuation and 
twined quasars) leads our analysis in the same direction 
of the raw theoretical model, that is, we get a Poisson-like 
distribution for all quasars. On the other hand, a finest 
analysis shows two peaks in the quasars distribution: the 
highest is due to non- twined quasars, whose light is not 
trapped and guided in the channels and the shortest is due 
to twined quasars . This fraction is in the range 7%-35%. 



8 



S. Capozziello & G. lovane: Quasar luminosity and twin efTects induced by 



filamentary and planar structures 







Fig. 6. Some examples of simulated 
refraction channels. 
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Fig. 7. Quasar redshift Z after at- 
tenuation process. 



Input number of quasars 


number of channels 


number of quasars 


Twined quasars 


Twined Q./Q. after atten. 


10" 


60 


1015 


109 


7% 


10" 


143 


977 


141 


14% 


10" 


484 


995 


260 


26% 


10" 


840 


1023 


337 


33% 


10* 


1154 


981 


461 


35% 



Table 1. Tab. 1. Simulation summary. 



6. Discussion and conclusion 

In this paper, we propose a mechanism by which it is pos- 
sible to explain two of major puzzles of quasar theory: the 
huge luminosity and the twin (or brother) objects. Both 
of them, till now, have invoked exotic explanations as the 
presence of a supermassive black hole inside the core ra- 



dius or lensing effects, where the position and the size of 
intervening deflector (usually a galaxy) has to be in same 
sense "fine tuned". By our cosmological refraction chan- 
neling effects both the issues are addressed in a natural 
way. 

The huge luminosity at far red-shift is nothing else but 
the effect of the channel which " draws" the luminosity of 
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Fig. 8. Scatter plot x — y plane 
of quasars after attenuation and 
twined quasars. 
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Fig. 9. Scatter plot x — z plane 
of quasars after attenuation and 
twined quasars. 



the source closer to the observer. In this sense, quasars 
could be just primordial (in same sense "ordinary") galax- 
ies whose luminosity in not dispersed as l/R"^ This fact 
means that the source brightness will turn out to be much 
stronger than the brightness of analogous objects located 
at the same distance (i.e. at the same redshift Z) and the 
apparent energy released by the source seems anomalously 
large. Besides, the large angular distance separation of 
twins and brothers is due to the geometry of (filamentary 
or planar) channels without posing "ad hoc" lens galax- 
ies between the source and the observer. This hypothesis, 
as shown, naturally explain the double peak luminosity 
function of quasars. 

However, what we have discussed here is nothing else 
but a toy model. In fact, the discrete nature of matter dis- 
tribution within the channels affect the light propagation 



(we have to take into account scattering points, dissipa- 
tion, absorbtion, and dispersion in a more realistic model). 

Depending on the distribution of discrete matter, we 
could have several places where the light could lake out. 

Furthermore, the distribution of our filamentary and 
planar structures has to be compared to the microwave 
background in order to preserve the overall homogeneity 
and isotropy of observed very large scale structure. 

All these issues deserve further studies and the com- 
parison to the observational data. 
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